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^ime is a great teacher, but untortunately it kills all its pupils^ Hector Berlioz = 
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Time plays an essential role in most real world fluid mechanics problems, 
e.g. turbulence, combustion, acoustic noise, moving geometries, blast 
waves, etc. 

Time dependent calculations now dominate the computational landscape at 
the various NASA Research Centers but the accuracy of these 
computations is often not well understood. 



Helicopter and Launch Vehicle Combustion and 
Tilt-Rotor Aerodynamics Analysis turbulence 
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In this presentation, we investigate error representation (and error control) 
for time-periodic problems as a prelude to the investigation of feasibility of 
error control for stationary statistics and space-time averages. 

o These statistics and averages (e.g. time-averaged lift and drag forces) 
are often the output quantities sought by engineers. 

o For systems such as the Navier-Stokes equations, pointwise error 
estimates deteriorate rapidly which increasing Reynolds number while 
statistics and averages may remain well behaved. 
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im^sa Motivating Example #1 : Cylinder Flow 
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Cylinder flow at Mach = 0.10, logarithm of |vorticity| contours 


Introduction 

Cylinder 

Row 

Computabilit 
of Outputs 

Nonlinear 

Conserva- 

tion 

Laws 

Space-time 

Prisms 

Space-time 

DG 

Error Rep- 
resentation 

Scalar 

transport 

Navier- 

Stokes 

Formulation 

Example 

Dual 

Problems 

Pariodic 

Cylinder 



Re=1000 



Re=3900 



Re= 10000 Re=50000 



o Quartic space-time elements t f * — 1#** ^ 
o 25K element mesh 

o Viscous walls only imposed on cylinder surface 
o Reynolds number based on cylinder diameter 
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Example: Backward facing step (Re=2000) 



Suppose J(u) is the streamwise velocity component averaged in cube in 
space and over a unit time interval, i.e. 

J(u) = [ f Uidx 3 dt 

•> 9 Jdxdxd 
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Hoffman and Johnson (2002) have computed solutions of the backward 
facing step problem using a FEM method with linear elements for 
incompressible flow. 

In velocity and pressure variables, ( V , p), the following error estimate for 
functionals is readily obtained in terms of the dual solution (^, <t>) 

\J(V,p)-J(V h ,p h )\ < C\\M\\Atro(V h ,Ph)\\ 

+ C\\[?mt?ro(V h ,p h )\\ 

+ C\m\Atr,(V h ,p h )\\ 

+ C||D^||||/7r,(Vh,p h )|| 


where r, are element residuals. 
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The following stability factors have been computed by Hoffman and Johnson 
(2002) for the backward facing step problem at Re=2000. 


d 

m 

livv'li 

l|v*|| 

11*11 

1/8 

124.0 

836.0 

138.4 

278.4 

1/4 

39.0 

533.4 

48.9 

46.0 

1/2 

10.5 

220.3 

16.1 

25.2 


These results clearly show the deterioration in computability as the box 
width is decreased. 


= ■OQ.O' 


im^sa Outline for the Remainder of the Talk 
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o Review the space-time discontinuous Galerkin (DG) FEM formulation, 
Reed and Hill (1973), LeSaint and Raviart (1974) and popularized for 
nonlinear conservation laws by Cockburn and Shu (1990). 

o Error representation and estimation for nonlinear hyperbolic systems 
with and without time 

o The space-time discontinuous Galerkin method for the compressible 
Navier-Stokes equations 

o Error representation and estimation for time periodic, and nearly time 
periodic Navier-Stokes cylinder flow 

o (Time Permitting) Recent work moving away from functional error 
representation/control towards T p -norm control. 


« □ ► « £5P ► «-=►«$► ■= ■OP.O' 


im^sa Nonlinear Conservation Law Systems 


Space- 

Time 

DG 

Tim Barth 


Introduction 

Cylinder 

Row 

Computability 
of Outputs 

Nonlinear 

Conserva- 

tion 

Laws 

Space-time 

Prisms 

Space-time 

DG 

Error Rep- 
resentation 

Scalar 

transport 

Navier- 

Stokes 

Formulation 

Example 

Dual 

Problems 

Pariodic 

Cylinder 


Conservation law system in R dx1 

Uj + divf = 0, u,f/ e R m /= 


Convex entropy extension 

Uj -b div F < 0, U, F-, e R 


*0 0.0 
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Existence of a convex entropy-entropy flux pair {l/, F} implies that the 
change of variable u^v symmetrizes the original quasilinear system (Mock 
( 1980 )) 


u, v v,t + f/,v v, x/ = 0 (implied sum, / = 1 ...d) 

SPD SYMM 

so that for smooth solutions 


v • (ll,f -f divf) = Uj + div F = 0 . 

with the symmetrization variables (a.k.a. entropy variables) calculated from 
v r = U, u and v f, v = F, v . 
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im^sa The Discontinuous in Time Approximation Space 
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o Natural setting for the discontinuous Galerkin (DG) method for 
hyperbolic problems 

o Utilized in the space continuous Galerkin least-squares method 
(Hughes and Shakib, 1988) 

o Often used in the discretization of parabolic problems (Douglas and 
Dupont, 1976) 

o Requires solving the implicit slab equations-no one said it would be 
easy! 



Discontinuous timeslab *■ a - orY , . ^ . 

. . Space-time prism element 

intervals 




Space-Time Discontinuous Galerkin Formulation 
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Piecewise polynomial approximation space: 

v h = {v h I VhlKx/n € (v k (K X / n )) m } 

Find V/, e V h such that for all w/, e V h 

A /— 1 

B(Vh< W/7)dg = B n (^h^h)DG = 0 , 
n = 0 


e n (v, w) DG = f v f -(U (V) • W,( + f'(v) w , Xi )dx dt 

Ji Wt Jk 

+ [ J2 [ w(x_)h(v(x_),v(x + );n)cfedf 

JrferJaK 

+ J (w(£ +1 ) • U(v(f" +1 )) - w(f") • u(v(f" ))) dx 


o u the conservation variables, v the symmetrization variables 
o ha numerical flux function, h(v_,v+;n) = -h(v+, v_; -n), h(v.v; n) =f(v) n 
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rt^sa Nonlinear Stability of Space-Time DG Formulations 
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Theorem E: Global space-time entropy inequality (Cauchy IVP): 

[ t/(u*(fi)) dx < [ U(u(v h (x,d)))dx< f U(u(v h (x,£)))dx 

Jn Jn ./n 

whenever the numerical flux satisfies the system extension of Osher’s 
famous “E-flux” condition 


[v]£ ■(h(v.,v + ;n)-f(v(0)).n)<O , V0g[O,1] , v(0) = v_ + 9[v]± 


o Several flux functions satisfy this technical condition when recast in 
entropy variables, e.g. Lax-Friedrichs, HLLE, Roe with modifications, 
etc. 
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rt^sa Nonlinear Stability of Space-Time DG Formulations 
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Suppose U'V remains bounded in the sense 


Vz^ 0 


and Theorem E is satisfied for the Cauchy IVP, then following L 2 stability 
result is readily obtained 


Z_2 Stability: 

I|U(V„(-, £) - u*(f° )|| 4(n) < {c^C o y /z ||u(v„(-, t °_ )) - u*(£)|| Mn) . 


«=►<=► = ■oq.O' 


Given a system of PDEs with exact solution u e R m in a domain ft, the overall 
objective is to construct a locally adapted discretization with numerical solution Uh that 
achieves 


o Norm control [Babuska and Miller, 1984] 

||u - u h \\ < tolerance 

o Functional output control [Becker and Rannacher, 1997] 

|J(u) - J(u h )\ < tolerance , J( u) : R m *-> R 

Example functional outputs: 

o Time-averaged lift force, drag force, pitching moments 
o Average flux rates through specified surfaces 
o Weighted-average functionals of the form 


for some user-specified weighting U/(x, t) and nonlinear function N(u) 



rt^sa Error Representation: Linear Case 
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Assume •) bilinear and J() linear. 

Primal Numerical Problem: Find U/, e such that 

S(U/ 7 ,w) = F(w) VweVJ. 

Auxiliary Dual Problem: Find <t> e V B such that 

B( w. <t>) = J(w) V w e V B . 

J(u) - J(u h ) = J(u -u h ) (linearity of J) 

= B( u — u h, <t>) (dual problem) 

= B ( u - U/ 7 , <J> - 7t/ 7<I>) (Galerkin orthogonality) 

= B( u. <t> - 7 t/j4>) - B(Uh , <t> - 7r/,<t>) (linearity of £) 

= F(<t> - ^0) - B(u h , 4> - tt/, 4)) (primal problem) 


Final error representation formula: 

J(U) - J(U/,) = F(4> - 7T^<t>) - 8(U/,. <t> - TTfc^) 
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NASA Estimating cfc> — 7177$: 


Space- 

Time 

DG 

Tim Barth 

Introduction 

Cylinder 

Flow 

Computabilit 
of Outputs 

Nonlinear 

Conserva- 

tion 

Laws 

Space-time 

Prisms 

Space-time 

DG 

Error Rep- 
resentation 

Scalar 

transport 

Navier- 

Stokes 

Formulation 

Example 

Dual 

Problems 

Pariodic 

Cylinder 


Various techniques in use for estimating <t> - 717,0 : 

o Higher order solves [Becker and Rannacher, 1998],[B. and Larson, 
1999], [Suli and Houston, 2002], [Houston and Hartman, 2002] 

o Patch postprocessing techniques [Cockburn, Luskin, Shu, and S uli, 
2003] 

o Extrapolation from coarse grids 
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Mean-value linearized forms: 

B( u.v) = £(U/,.v) + tf(u - U/,, v) Vve V B 

J(u) = J(u h ) + J(u-u h ), 

Example: B(u, v) = ( L(u ), v) with L(u) differentiable 

r u B r u B dL 

= J ^ (u(9 )) dO (u B -u A ) = L,u ■ (u B - u A ) 

with u(0) = ua + (ub - ua) 9. 

H(u,w) = B(u h ,)N) + (E'U • (u - U/,),w) 

= B(U h ,\N) + B(\J -u h .\N) W g V B 
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rt^sa Error Representation: Nonlinear Case 
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Semilinear form B ( •, •) and nonlinear J(-). 

Primal numerical problem: Find e V B such that 

#(U/,,w) = F(w) VweV B . 

Linearized auxiliary dual problem: Find <t> e V B such that 
#(w, <t>) = J(w) VweV B . 

J(u) - J( Uh) = J( u - u h) (mean value J) 

= B(u - U/ 7 , <t>) (dual problem) 

= B ( u — U/ 7 . — ^0) (Galerkin orthogonality) 

= #(u, <t> — n h <$>) — B{u h , <t> — 7T/,<1>) (mean value #) 

= F(<J> - 7r/,<t>) - #(11/,, <t> - 7r/,<t>), (primal problem) 


Final error representation formula: 

J(U) - J(U h ) = F(4> - n h <t>) - B(U h . 4) - n h <t>) 




Space-time error representation formula 

A /-1 

Edg(V/7. w) - Fdg( 4> - nh<t>) = YY ^DG,O n ( v /7. $ - 7r /7 <t> ) - ^DG.Q"^ ~ * h &) 

n= 0 Q" 


Stopping Criteria: 


\J(u) - J(u h )\ = 


N-^ 

Y Y ^DG.Q"( v / 7< - ^DG ,O n ( < * > - ^Ti^) 

n=0 0° 


Refinement/Coarsening Indicator: 


A /-1 


\J(U) - J(U/,)| < ^ Y \ B DG,Q n (^h^ - Kh^) - F D G,O n {^ - TT/J^) 
n = 0 Q n " ' 


refinement indicator 


o This provides a unified framework for both stationary and time dependent 
problems 


Example Error Representation in Space-Time 
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From the error representation formula, weighted estimates are obtained in space-time 

N 

J(U) - J(U h ) = X! (( r,7? ^ “ n h^)Qn + (j/7. $ - 7 ThtyaQn) 
n= 0 Q n 

where denotes the residual on element interiors 

r/, = u M + div(f(u/,)) . 

and denotes one of four possible jump terms 

'f(n;u/ 7 (x_))-h(n;u/ ? (x_),u /? (x + )), dQ n \ r, t^O 
_ I f(/?;U/ J (x_))-h(/?;U/ 7 (x_),g(x + )), dQ^nr 

to = 1 (u h (x.t + )-u h (x,t-)), ao"n[r]l 

(u„(x,f)-u 0 (x)), d&>, t = 0 
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im^sa Example: A Scalar Time-Dependent PDE 
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Circular transport, A = (y, -x), of bump data 


u t + A • Vi/ = 0 , 

u(x, 0) = v|/(1/10 ; x - Xq) , 


x e [— 1 , 1 ] 2 
X 0 = (7/10.0,0) 




Primal numerical problem Convergence, || u - ^/ 7 ||/. 2 (n x [o, 7^) 


«=►<=► = ■OQ.O' 


f ' - 

ti%$t \ Space-Time DG Method 


Space- 

Time 

DG 

Tim Barth 


Introduction 

Cylinder 

Flow 

Computability 
of Outputs 

Nonlinear 

Conserva- 

tion 

Laws 

Space-time 

Prisms 

Space -time 
DG 

Error Rep- 
resentation 

Scalar 

transport 

Navier- 

Stokes 

Formulation 

Example 

Dual 

Problems 


Example: Circular transport of bump data, A = (y, - x ) 

l/| + A • Vu = 0 , X€[-1.lf 

3K element mesh 





in space-time 


V 2 in space -time 


Periodic 
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im^sa Example: A Scalar Time-Dependent PDE 
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A functional is chosen that averages the solution data in the space-time ball of radius 
1/1 0 located at x c = (1 /2, 1 /2, 1 .05) in space-time 

J(u)= f f U/(1 /1 0; x - x c )udxdt 
Jo Jn 

o 

J ( u ) - J (Uh) = F DG,Q n {^ ~ ~ %G,0"( v /»> ^ 

n=N—'\ K 
0 

M(U) — J(u h )| < Z:i f DG.o n ( <t>_7rh<t> ) _e DG.O' , ( v h <t) - 7r h < t>)| 

n=N—1 K 




Software Implementation and extension to the Navier-Stokes 
Eqns 
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Space-Time FEM: 

o Implements the discontinuous Galerkin discretization in entropy variables. 

o Unconditionally stable for all time step sizes 

o Parallel implementation using overlapping domain decomposition and ILU 
preconditioned GMRES subdomain solves. 

o Solves both the primal numerical problem and the jacobian linearized dual 
problem arising in space-time error estimation. 

o High-order accuracy demonstrated in both and CSB 

o DG extension to the compressible Navier-Stokes equations using the symmetric 
interior penalty method of Douglas and Dupont, 1976) as described in Hartmann 
and Houston (2006) 
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Space-Time DG Formulation for the Navier-Stokes Eqns 


Find \ih £ such that 


Space-time 

Prisms 


Space-time 

DG 
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Formulation 


with 


N - 1 

Bdg(V/>, W) = Y1 Bdg, 0"(V/„ W) = 0 , Vw g v“p 

n=0 K 

Sdg ,cy»( v - w ) = 


+ f in •(*- )■ Ml* *+.»->-<» F(”(»-)) 

+ f„ r L «<*- ) • "/ (if wa " - if )(»- )) dt dx 

+ /-nr_ Im "<*- > • (h(n: .»-)-'>/ f (*-»<* * 

- /«rnr„ /* W(X - ) ■ n/ ,F?iS(flw) " F ' iS(v - ))d,ds 
-L\rfinl* x - )ln rf' ]X *- d,dx 

+ I S K\rfl n 

+ / 0/Cn r D - v (*-» • niMij(x_)vi fXj (x_)dtdx 

- J dKV f in <«f/ h >xl w ( x - » • n i n i M ij yfx + _ <* * 

- feKnr D Jin {Kf ^ /h)x - • °/ n j M ij < g D - v(x_ ))) dtdx 
+ J K w ( f +) • t u ( v )l[i <* + ^ ^ 0 W (^) • (W-)) - Uo)d>r 


Pariodic 

Cylinder 


im^sa Primal-Dual Problems in Fluid Mechanics 
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Subsonic Euler flow 

M = .10, 5° AOA 
Primal Mach contours 
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Transonic Euler flow 
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M = .85, 2° AOA 
Primal density contours 




Lift force functional 

Dual x-momentum 
contours 



Lift force functional 

Dual density contours 
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Viscous cylinder flow 

M = .15, Re = 300 
Primal vorticity contours 




Drag force functional 

Dual x-momentum 
contours 


«=» <=► = ■oq.O' 


Periodic 

Cylinder 


im^sa An Application of Error Estimation and Adaptive Error Control 
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Example: Euler flow past multi-element airfoil geometry. M = .1 , 5° AOA. 


lift coefficient 
(error representation) 

lift coefficient 
(error control) 

refinement 

level 

# elements 

equivalent uniform 
refinement 
# elements 

5.1 56 dz .1 47 

5. 156 ±.346 

0 

5000 

5000 

5.275 ± .018 

5.275 ± .076 

1 

11000 

20000 

5.287 ± .006 

5.287 ± .024 

2 

18000 

80000 

5.291 ± .002 

5.291 ± .007 

3 

27000 

320000 



Error reduction during mesh adaptivity Adapted mesh (18000 elements) 


Periodic 

Cylinder 
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Computing dual (backwards in time) problems looks expensive in terms of both 
storage and computation 


Tremendous simplification arising for periodic flow problems with period P when 
phase-independent functionals are utilized, e.g. mean drag 
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Computing dual (backwards in time) problems looks expensive in terms of both 
storage and computation 

o Storage of the primal time slices for use in the locally linearized dual problem. 


Tremendous simplification arising for periodic flow problems with period P when 
phase-independent functionals are utilized, e.g. mean drag 
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Computing dual (backwards in time) problems looks expensive in terms of both 
storage and computation 

o Storage of the primal time slices for use in the locally linearized dual problem. 

o Approximation of the infinite-dimensional dual problem for the backwards in time 
dual problem. 


Tremendous simplification arising for periodic flow problems with period P when 
phase-independent functionals are utilized, e.g. mean drag 
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Computing dual (backwards in time) problems looks expensive in terms of both 
storage and computation 

o Storage of the primal time slices for use in the locally linearized dual problem. 

o Approximation of the infinite-dimensional dual problem for the backwards in time 
dual problem. 


Tremendous simplification arising for periodic flow problems with period P when 
phase-independent functionals are utilized, e.g. mean drag 

o Functional independent of the startup transient 
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Computing dual (backwards in time) problems looks expensive in terms of both 
storage and computation 

o Storage of the primal time slices for use in the locally linearized dual problem. 

o Approximation of the infinite-dimensional dual problem for the backwards in time 
dual problem. 


Tremendous simplification arising for periodic flow problems with period P when 
phase-independent functionals are utilized, e.g. mean drag 

o Functional independent of the startup transient 

o Only a small number of periods of the primal problem need be stored or 
recreated. 
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Cylinder flow at Mach = 0.10, logarithm of |vorticity| contours 



Re=300 Re=1000 

Task: Represent and estimate the error in the mean drag force coefficient 
o Solve the primal problem using linear space-time elements 
o Construct a smooth phase invariant functional measuring the mean 
drag force coefficient 

o Solve the dual (backwards in time) problem using quadratic space-time 
elements 

o Calculated the estimated functional error and compare with a reference 
calculation using cubic elements 
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J dra g (u) = / 0 r / rwt|| (Force • k«g) V(t) dx dt 
Example: Cylinder flow at Re=300 
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Dual problem at Re=300 
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rt^sa Mean Drag for Cylinder Flow at Re=1 000 
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Adapted mesh from element indicators averaged over a period P 




Coarse mesh (12K elements) 2 level refined mesh (20K elements) 


« □ ► « £5P ► «-=►«$► ■= 


im^sa Non-Periodic Cylinder Flow 


Space- 

Time 

DG 

Tim Barth 

Introduction 

Cylinder 

Row 

Computabilit 
of Outputs 

Nonlinear 

Conserva- 

tion 

Laws 

Space-time 

Prisms 

Space-time 

DG 

Error Rep- 
resentation 

Scalar 

transport 

Navier- 

Stokes 

Formulation 

Example 

Dual 

Problems 

Periodic 

Cylinder 


Cylinder flow at Re=3900 and Re=10000. 

o Choosing measurement problems that are not genuinely stationary 
produces rapidly growing dual problems and dependency on the initial 
data. 



Dual solution corresponds to the average drag force over 3 approximate 
“periods”. 



a Concluding Technical Remarks 
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o Including time as “just another dimension” has many merits 
o Arbitrary order approximation 
o Provable non-linear stability 
© Simplified space-time error estimation 

o But it also comes at a price 

© Increased arithmetic operations 
© Increased memory storage 
© More complex code implementation 

o Error representation/estimation results presented today barely scratch 
the surface 

o Error control for general transient problems, 
o Dual problems in the presence of flow bifurcations 
o Computability and deterioration of functionals with increasing Reynolds 
number 

© Computer memory and storage constraints. 
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Schematic of Ringleb flow Iso-Density contours 



Discontinuous Galerkin 
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im^sa Example: A Scalar Time-Dependent PDE 
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Circular transport, A = (y, —x), of bump data 


Ut -|- A • V u — 0 , 

U{x, 0) = V(1/10;X — Xo) , 


x e [— 1 , 1] 2 

x 0 = (7/10.0.0) 




Primal numerical problem Convergence, || u — ^/7lk 2 (nx[o,r]) 


